Spin-dependent thermoelectric transport through double quantum dots 
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We study thermoelectric transport through double quantum dots system with spin-dependent in- 
terdot coupling and ferromagnetic electrodes by means of the non-equilibrium Green function in the 
linear response regime. It is found that the thermoelectric coefficients are strongly dependent on the 
splitting of interdot coupling, the relative magnetic configurations and the spin polarization of leads. 
In particular, the thermoelectric efficiency can achieve considerable value in parallel configuration 
when the effective interdot coupling and tunnel coupling between QDs and the leads for spin-down 
electrons are small. Moreover, the thermoelectric efficiency increases with the intradot Coulomb 
interactions increasing and can reach very high value at an appropriate temperature. In the pres- 
ence of the magnetic field, the spin accumulation in leads strongly suppresses the thermoelectric 
efficiency and a pure spin thermopower can be obtained. 
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I. INTRODUCTION 

In recent years, quantum dot (QD) systems, as 
a potential artificial thermoelectric material, have at- 
tracted much attention and their thermoelectric prop- 
erties have been widely studied both experimentallyii^— 
and theoretically^ ~^^'^^~ — . Due to deviation from the 
Wiedemann-Franz law and suppression of the phonon 
contributions, thermoelectric efficiency in QD systems 
can be significantly enhanced^i^i^. On the other hand, 
many interesting thermoelectric effects have been ob- 
served or predicted in single QD and double quantum 
dots (DQDs) systems in the Coulomb blockade (CB) 
regime^"— and Kondo regimei^"— . For example, the 
thermopower and the thermal conductance as a func- 
tion of gate voltage present CB oscillations''. The tran- 
sient thermopower of a single QD system can be strongly 
enhanced and modified under a time-dependent gate 
voltage^. In a single QD with multilevelsi^ or a serial- 
coupled DQDsi^ii^, the thermal conductance can be 
strongly enhanced at high temperature due to the bipo- 
lar effects, leading to strong suppression of thermoelec- 
tric efficiency. When interference effects are considered 
in parallel coupled DQDs system, the thermoelectric ef- 
ficiency can be significantly enhanced by Fano effect and 
Coulomb interactioniir— In particular, the strong de- 
pendence of thermopower on temperature and gate volt- 
age provide a sensitive way to detect and understand the 
Kondo correlations2ii2i. 

The recent experimental observation of the spin See- 
beck effect^"— inspires the research on spin effects in 
energy transport in QD systems^iSii^Sri^ and opens a new 
view of thermal manipulation of spintronic devices and 
the design of spin current generator. Spin-dependent 
thermoelectric effects in the single QD and DQDs system 
connected to two ferromagnetic leads with coUinear or 
noncoUinear magnetic configurations is widely studied in 
the linear and nonlinear response regim o^'^'^'^i'^°~ '— . The 
results indicate that the thermoelectric coefficients are 



strongly dependent on the spin polarization, the relative 
magnetic configurations of ferromagnetic leads and the 
asymmetry of dot-lead coupling. As mentioned by Jozef 
Barnas^i^i^, due to long spin-relaxation time, the spin 
accumulation in leads induces spin-dependent splitting of 
the chemical potential in the leads. As a result, the spin 
efficiency can be larger than the charge efficiency and a 
pure spin thermopower can be obtained^*^. Up to now, 
most earlier theoretical works about spin effects in heat 
transport mainly focused on single QD or parallel cou- 
pled DQDs system. More recently, people start to study 
thermoelectric effects in a serial-coupled DQDs system 
coupled to ferromagnetic leads^^. Even so. Still much 
more work is needed to understand interplay between 
spin effect and heat transport in serial-coupled DQDs 
system, particularly in the presence of the magnetic field 
or the spin accumulation in leads. 

In the present work, we consider a lateral DQDs with 
spin-dependent interdot coupling connected to two fer- 
romagnetic electrodes. Spin-dependent interdot coupling 
can be achieved by applying a static magnetic field on the 
tunneling junction between two QDs and induces the lev- 
els spin-dependent splitting. The tunable system prop- 
erties further enrich the thermoelectronic and thermo- 
spin transport properties in DQDs system. By using 
the nonequilibrium Green's function technique, thermo- 
electric and thermo-spin coefficients are obtained. It is 
found that in parallel configuration the thermoelectric 
efficiency can achieve considerable values around spin- 
down resonance levels when the effective interdot cou- 
pling and tunnel coupling between QDs and the leads for 
spin-down electrons are small. On the other hand, in the 
presence of the magnetic field, the spin accumulation in 
leads strongly suppresses the thermoelectric efficiency. In 
particular, the thermoelectric and thermo-spin efficiency 
are strongly enhanced by the intradot Coulomb interac- 
tions and can reach very high value at an appropriate 
temperature. Moreover, a pure spin thermopower can 
also be obtained. 
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The rest of this work is organized as foUows. In Sec. II, 
we introduce the model of lateral DQDs system with 
spin-dependent interdot coupling and derive the basic an- 
alytical formulas. In Sec. Ill, we present the correspond- 
ing numerical results about the influence of the related 
parameters on the thermoelectric and thcrmospin prop- 
erties in the present device. Finally, In Sec. IV we sum- 
marize the work. 



II. MODEL AND FORMALISM 

We consider a system of lateral DQDs with spin- 
dependent interdot coupling and two ferromagnetic leads. 
The magnetizations of the leads are assumed to be 
collinear, and can be either in the parallel or antiparal- 
lel configuration. The Hamiltonian describing the sys- 
tem can be decomposed into three parts, i.e., H = 
Hd + Ha + Ht- The first term describes the isolated 
DQDs and is given by. 



we have 
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with = T1^^ and TJ^ = r^2fl,- When the spin po- 
larization of the leads is considered, the matrix-element 
may be assumed in the form: F^^ = Fl(1 + So-Pl) and 
Fj^ — Ffl(l ± SaPn) with the sign -|- (— ) for parallel 
(antiparallel) alignment of the magnetization in two elec- 
trodes. Pa is related to the lead's spin polarization. F" 
describes the tunneling coupling between the DQD and 
the a-th electrode for pa = 0. In the following calculation 
we assume Pl = PR = P, and F^ = F^,; — F. 

By using nonequilibrium Green's function technique, 
the electronic and heat current from the a (L, R) lead 
flowing into the QDs can be written in the formsi^i^i^S: 



/^H]r,H (3) 



W - Aia) [fa (w) - //3 (w)] T„ (w) (4) 



where the operator dl^^ (dma) creates (annihilates) an 
electron with energy Ema and spin a in QD m. We as- 
sume that the level spacing is very large so that we can 
only consider a single energy level in each QD. U is in- 
tradot Coulomb interaction constant and Uma = d\^^d.ma 
is the particle number operator. The last term sum in 
Eq.(l) describes the spin-dependent tunneling between 
two QDs, where 5a is defined as: 5a ~ ^ for a =t 
or 5a — —1 for a =\r. This spin-dependent tunneling 
strength can be induced by a static magnetic field ap- 
plied to the barrier region between the two QDs. When 
electron passes through the potential barrier region, the 
electron with spin parallel to the magnetic field have a 
larger tunneling probability than one with spin antipar- 
allel to the magnetic field due to Larmor precession of 
the electron spin inside the potential barrie r^^i'^^ . 

The second term Ha of Hamiltonian describes the non- 
interacting electrons in the leads and is written as Ha — 
T,a=L,R T,ka ^kaAaaCkaa , whcrc 4„^(cfcc«^) Creates (an- 
nihilates) an conduction electron with energy Ska, rao- 
mentum k and spin a in the a electrode. The third term 
of Hamiltonian describes the coupling between the DQDs 
and the electrodes and reads Ht — J^kai^Liacli^a^^cr + 
yR2iyc\j^^d2a +H.c), where Vama IS the hopping matrix- 
element between the QD m and the electrode a, which 
is assumed to be independent of momentum k. In the 
following calculation, we define the line-width matrix 
as T^nrna = ^Q«cr Sfc 27r(5(a; - £fca ) (a = i,i?,and 
Vl2(7 = Villa) and in a wide-band limit, FJ^^^. is inde- 
pendent of the energy. Using the matrix representation. 



where f3 ^ a, fa = [e("-'^°)/'=«^" + is the Fermi 
distribution function for the a electrode with chemical 
potential /Iq, temperature Ta and Boltzmann constant 
ks- Ta (w) is the transmission coefficient corresponding 
to spin channel a and is given by 



(5) 



Here GJJ.(w) and G° (w) are, respectively, the retarded and 
advanced Green's function in the DQD in uj space, with 
G°(a;) — [G^{uj)]'' . Employing the equation of motion 
and adopting the Hartree-Fock truncating approximation 
for higher order Green's functions, the retarded Green's 
function can be obtained as from the Dyson equations 



G22A^) 



(6) 

Where SJ^ = —i [Tp + FJ] is the retarded self-energy in 
the wide-band limit and gai^) is the retarded Green's 
function for the isolated DQD (without coupling to elec- 
trodes) which is expressed as 



Cia t + 5a At 
t + 5aAt C2a 



(7) 



where Cma (w - Sma - U){lO - eraa)l{^ - £mcr - U + 

Urima) and Uma, the averaged occupation number of 
electrons with spin a in the m-th quantum dot, can be 
calculated self-consistently by using the formula Uma = 

- -ij i^G<a,msi^), where G<(w) is the 
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lesser Green function in lo space and can be obtained by 
the Keldysh equations G<{uj) = Gl{uj)Y.< {uj)G%{uj) with 
the lesser self-energy E<(a;) = i[fL{w)Tl + fRiu;)^^^]. It 
should be pointed that the Hartree-Fock decoupling ap- 
proximation scheme is reasonable when the temperature 
during our calculation is higher than the Kondo temper- 
ature. 

In the linear response regime, the spin accumulation in 
electrodes can be neglected when spin relaxation time is 
sufhciently shor t^'^^i'^^ . In this case, the voltage induced 
by temperature gradient can be treated as independent of 
spin. The electric and heat currents through the system 
can be written in the form 



/ ^ e^LoAV 
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with L^^-E^yi^iuj-f^r m^T. (c.). In the 
following calculation we assume the chemical potential 
Ml = I^Ri and temperature = Tr, thus the Fermi 
distribution function fai^jj) — fi'-^)- The electronic and 
thermal conductances can be calculated by the formula 

G = c^Lq and k = ^[^2 — j^], respectively. The ther- 
mopower is defined as thermoelectric voltage AT^ induced 
by a temperature gradient AT under the condition of 
/ = 0, i.e., S = ^ = --^T-^- The efficiency of heat- 
electricity conversion of the system is described by the 
dimensionless figure of merit (FOM) ZT = ^"^^^ , where 
Kph denotes the thermal conductance corresponding to 
the phonon contribution, which is suppressed in such 
DQD-electrodes structure, and is ignored in the present 
paper. 

If spin relaxation time of systems is sufficiently 
long, the spin accumulation in electrodes should be 
considere d^'^°i'^^ , which leads to that the chemical po- 
tential difference AV^ between the two leads is related 
to the spin channel. Thus, in this system one can de- 
fine a spin bias AYs induced by temperature gradient 
via the formula AV^ — AY -f 8a AY^. Corresponding 
spin-resolved electric and heat currents can be expressed 
as 



,^1 



e'LlAYa^e^AT, 



E 



eLlAYa + ^AT 



(10) 
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with LI = -ll^^iu,-^,f^-M^Ta{^) in ^ 
0,1,2). The charge and spin conductance are calcu- 
lated as G = e^\L\ ^ L;^] and Gs = e^\L\ - L;^], re- 
spectively, while thermal conductance is given by k = 
— ^[-^2 27^] 5 which is identical to the value with- 
out spin accumulation in such a system. One may in- 
troduce from Eq.(lO) the spin-dependent thermopower 



Sc — ^f- — —-^jy- in spin channel a under the con- 
dition of current I„ vanishing^. Then the charge and 
spin thermopower are defined as Sc = = ^ [S^^ -\- S\\ 
and Ss — = \ [S^ — Si\ , respectively. Corresponding 
charge and spin figure of merit are defined as ZTc — 
and ZT^ — 



GSiT 



respectively, which describe heat- 
to-charge- voltage and heat-to-spin- voltage conversion ef- 
ficiency of the present system when spin accumulation is 
considered. ZTg has the form of absolute value because 
the spin conductance may be negative in some regions. 



III. NUMERICAL RESULTS AND DISCUSSION 

We now numerically study spin-dependent thermoelec- 
tric effects in a DQDs system based on the above formula, 
which will be presented according to three premises: 
without spin accumulation in the leads and U — 0] with- 
out spin accumulation in the leads and C/ 7^ 0; with spin 
accumulation in the leads. For the sake of simplicity, we 
assume two dots have same bare energy levels Sma and 
are spin degenerate, i.e. Sma = £, which can be tuned by 
gate voltage in experiment. The relevant parameters are 
chosen as: = ^ir = 0, Tl ^Tr = T,T = 0.1 meV, 
t = 2 meV. 



without spin accumulation and [7 = 



Consider the case that the magnetizations of the leads 
are in parallel configuration, there is no spin accumula- 
tion in the leads and the intradot Coulomb interactions 
are ignored. The tunnel coupling between two QDs is 
spin-dependent due to the external magnetic field, such 
a coupling DQDs can be equivalent to a single QD with 
four spin-dependent levels, Ei^f — e + t + At, = 
e-t- At, E31 ^ e + t-At, and E41 ^ e-t + At, which 
corresponds to poles of the Green function .gj('^) given 
by Eq.(7). These levels are identified as Euf (u=l,2) oc- 
cupied by spin-up electrons and Ed^ (d=3,4) occupied by 
spin-down electrons, and satisfy Ei^ > E^i > E41 > E2-\ 
for chosen parameters. 

The influence of splitting At of interdot coupling on the 
thermoelectric transport is shown in Fig. 1. Four reso- 
nance peaks appear in the electrical conductance spec- 
trum when the spin-dependent level crosses the Fermi 
level [fjL = 0), see Fig. 1(a). It can be clearly seen 
that the variety of At only changes the electrical con- 
ductance's resonance positions and do not influence its 
magnitude. It is found that this hold true for the ther- 
mal conductance n at low temperature. However, the 
behavior of k is more complicated at relatively high tem- 
perature (shown in Fig. 1(b)), and the thermopower S 
and the figure of merit ZT are more sensitive to the vari- 
ety of At. The antisymmetry, oscillation and sign rever- 
sal of S also are obtained from Fig. 1 (c) . As shown in Fig. 
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1(d), ZT presents typical double-peak structure. When 
Ai increases, for the two peaks of ZT located around 
e = Ed^, not only their magnitudes augment but also 
their degree of asymmetry increases. In particular, ZT 
exhibits a behavior similar to crossover and can reach 
very high value when Ai is close to t (see the inset of 
Fig. 1(d)). This is because when At is close to t, the 
levels E'aj. and E/^\^ get much closer to each other and the 
spin-down resonance channels are significantly blockaded 
due to the strong reduce of the interdot tunneling for 
spin-down electrons. While with At increasing the side 
peaks located at e = iS^t only shift their positions with 
unchanged the maximum. 

Next, we study the influence of spin polarization p on 
the thermoelectric transport. The relevant thermoelec- 
tric quantities versus dot's level position e for different 
spin polarization p are presented in Fig. 2. When p 
increases, the electrical and thermal conductance peaks 
(the thermopower peaks) located at e = Eu-\ increase 
(decreases) while the electrical and thermal ones (the 
thermopower ones) located at e = Ed\. decrease (in- 
creases) because Pj^, increases and F^^, Fj^ de- 
creases. These properties of G, k and S determine the 
behavior of ZT . As shown in Fig. 2(d), with increasing 
p the ZT near e = Eu-\ is suppressed while the ZT near 
e — Edi is intensively enhanced. The dependence of ZT 
on spin polarization in the leads and on spin-dependent 
levels may be useful for designing spin-dependent ther- 
moelectric devices. However, it should be noted that 
when the spin polarization in leads is very high, for ex- 
ample, p > 0.95, the tunneling rate for spin-down elec- 
trons is close to zero leading to significant suppression 
of transport through the spin-down levels and thus the 
corresponding ZT peaks rapidly decrease with increase 
of p, which can be clearly seen from the inset (2) of Fig. 
2(d). 

When the magnetic configuration changes from parallel 
to antiparallel, the electrical conductances are suppressed 
for four resonance channels, but the changes of the ther- 
mal conductance k and the thermopower S are dependent 
of the resonance levels (see Fig. 3). k {S) lowers (height- 
ens) for the spin-up resonance channels and heightens 
(lowers) for the spin-down resonance channels. As a 
result, ZT in antiparallel configuration presents higher 
peaks around the spin-up resonance channels and lower 
peaks around spin-down resonance channels comparing 
with in parallel configuration. For antiparallel config- 
uration, the tunneling rates between the leads and the 
QDs satisfy rJ^F2^ = F^^Fj^, which leads to the trans- 
mission coefficient (w) = \Gai2{'^)\'^^'iL^2R spin- 
independent. Then the thermoelectric quantities have 
the same intensities at different resonance points. Al- 
though the electric and thermal conductance experience 
an decreasing trend with increasing p, the thermopower 
and the figure of merit are independent of p. Based on 
the discussions above, in order to obtain high thermoelec- 
tric efficiency one should tune the resonances points to 
spin-down levels and take ^""0.9 in parallel configuration 



for fixed temperature and splitting of interdot coupling. 
B. without spin accumulation and U 



Now we focus on the case taking into account the in- 
tradot Coulomb interactions. Under the Hartree-Fock 
approximation, such a coupling DQDs system can be 
equivalent to a single QD with eight energy levels, 
in which four levels are occupied by spin-up electrons 
and other four are occupied by spin-down electrons^. 
To determine the eight spin-dependent effective en- 
ergy levels one need self-consistently solve the equation 
|[g^(w)]"i| = due to J7 7^ 0. Fig. 4 presents the 
relevant thermoelectric quantities versus the dots level 
position e with different temperature and U = 5. At 
low temperature eight resonance peaks of G and k can 
be clearly seen from Fig. 4(a) and (b). The peaks 1, 
4, 5, 8 (2, 3, 6, 7) correspond to the electron transport 
through spin-up (down) resonance channels. With the 
temperature increasing the conductance peaks become 
wider and lower due to broadening of the Fermi distribu- 
tion function. More interestingly, at high temperature, 
the thermal conductance presents five high peaks due to 
the bipolar eSecta^^^. 

S and ZT experience a nonmonotonous variation trend 
with increasing temperature (sec Fig. 1(c), (d)). In par- 
ticular, at high temperature, two high peaks of ZT ap- 
pear near symmetry point e = —U/2, resulting from a 
suppression of k and an enhancement of S. Comparing 
with the case in the absence U, it can be found that 
Coulomb interactions shift the positions of the peaks 
of ZT and strongly enhance ZT. As shown in the in- 
sert of Fig. 4(d), The maximum of ZT increases with 
Coulomb interactions increasing and can reach very high 
values above 30 at an appropriate temperature, which is 
very useful for the thermoelectric energy conversion de- 
vices. It should be mentioned that with varying the rela- 
tive magnetic configurations and the spin polarization of 
leads, the relevant thermoelectric quantities exhibit same 
variation trend as that without Coulomb interactions ex- 
cept for peaks around e = —U/2, which are unaffected. 

C. with spin accumulation 

For the systems with long spin relaxation time, the spin 
accumulation in electrodes should be considered^'^'^. 
Fig. 5 presents the thermoelectric coefficients versus e in 
the presence of spin accumulation. It is found by com- 
paring Fig. (3) with Fig. (5) that the spin accumulation in 
electrodes has no influences on the electrical conductance 
and the thermal conductivity (not shown) but strongly 
influences on the thermopower and the flgure of merit. 
The peaks of Sc and ZTc are much lower than the ap- 
propriate peaks in S and ZT. These properties are quite 
different from those in the absence of external magnetic 
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field, i.e. tlie splitting of interdot coupling At = 0. Since 
the levels are strongly spin-dependent in such system, the 
spin FOM experiences the same variation trend as charge 
FOM with varying the relative magnetic configurations 
and the spin polarization of leads, and the behaviors are 
similar to the case without spin accumulation discussed 
above. 

Some novel properties also appear, for instance, in 
the absence of Coulomb interactions, Gg = — G at reso- 
nance points e = E^^ while Gg = G at resonance points 
e = i?dt (see Fig. 5(a)). More interestingly, for certain 
£ (pointed by arrow and being around e = ±t in Fig. 
5(b)), the charge thermopower vanishes while the spin 
thermopower can be finite because the thermopowers for 
different spin have same value but opposite sign, which 
means that a pure spin current without charge current 
may be produced by a temperature gradient in our sys- 
tem. This feature can also be obtained in the spectrum 
of figure of merit in Fig. 5(c). When the Coulomb in- 
teraction is considered, near the symmetry point ZTc ex- 
hibits two very high peaks while ZTs completely vanishes. 
This is because the thermopower for different spin have 
the same value and sign leading to the enhancement of 
charge thermopower and the suppression of the spin ther- 
mopower. Even so, ZTc and ZTs are strongly enhanced 
by the Coulomb interactions, and ZT^ can be larger than 
ZTc for selected dots level positions in Coulomb block- 
ade regime. This implies our proposed system is a good 
heat-to-spin-voltage converter. 



IV. CONCLUSION 

In summary, we studied spin effects in thermoelectric 
transport through a lateral DQDs system with spin de- 
pendent interdot coupling and ferromagnetic electrodes. 
It is found that in parallel configuration the thermo- 
electric efficiency can achieve considerable values around 
spin-down resonance levels when the effective interdot 
coupling and tunnel coupling between QDs and the leads 
for spin-down electrons are small. On the other hand, 
in the presence of the magnetic field, the spin accumu- 
lation in leads strongly suppresses the thermoelectric ef- 
ficiency. The thermoelectric and thermo-spin efficiency 
are strongly enhanced by the intradot Coulomb interac- 
tions and can reach very high value at an appropriate 
temperature. Moreover, In such a DQDs system a pure 
spin thermopower can be obtained. 
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Figure Caption: Fig. 1 (Color online) Thermoelectric 
coefficients (a) G, (b) k, (c) S, and (d) ZT as a function 
of dot's level position e for varied splitting At in parallel 



configuration with ksT = 0.1 mcV. At is measured in 
meV. Other parameters are taken as: t = 2 meV, U = 
meV, F = 0.1 mcV, and p = 0.5. The inset in (d) depicts 
ZT as a function of e and At. 

Fig. 2 (Color online) Thermoelectric coefficients (a) G, 
(b) K, (c) S, and (d) ZT versus the dot's level position e 
for different spin polarization p in parallel configuration 
with ksT = 0.01 meV. At = 1 meV, other parameters 
arc the same as in Fig.l. The inset (1) in (d) is an en- 
largement of the region near £ = — 1 meV. The inset (2) 
depicts ZT as a function of spin polarization p with fixed 
e = -1.03 mcV. 

Fig. 3 (Color online) Thermoelectric coefficients (a) G, 
(b) K, (c) S, and (d) ZT versus the dot's level position e 
in different configuration: antiparallel configuration and 
parallel configuration with p = 0.5. Other parameters 
are the same as in Fig. 2. Inset: enlargement of the 
region near e = — 1 meV. 

Fig. 4 (Color online) Thermoelectric coefficients (a) G, 
(b) K, (c) 5*, and (d) ZT versus the dot's level position 
e for varied temperature T in parallel configuration with 
U = 5 meV and p = 0.5. The energy parameter fcsT is 
measured in meV. Other parameters are taken as in Fig. 
2. 

Fig. 5 (Color online) (a) Electrical conductance G and 
spin conductance Gg, (b) thermopower S and spin ther- 
mopower Sg , (c) (d) charge figure of merit ZTc and spin 
figure of merit ZTs versus e in parallel configuration with 
p = 0.5. The Coulomb interactions U = meV in (a) 
(b) (c) and U = 5 meV in (d), Other parameters are the 
same as in Fig. 2. 
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